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Abstract 

Discrete-time models are the traditional approach for capturing population dynamics of a host- 
parasitoid system. Recent work has introduced a semi-discrete framework for obtaining model update 
functions that connect host-parasitoid population levels from year-to-year. In particular, this framework 
uses differential equations to describe the hosts-parasitoid interaction during the time of year where they 
come in contact, allowing specific behaviors to be mechanistically incorporated into the model. We use 
the semi-discrete approach to study the effects of host-feeding, which occurs when a parasitoid consumes 
a potential host larva without ovipositing. Our results show that host-feeding by itself cannot stabilize 
the system, and both the host and parasitoid populations exhibit diverging oscillations similar to the 
Nicholson-Bailey model. However, when combined with other stabilizing mechanisms such as density- 
dependent host mortality or density-dependent parasitoid attack rate, host-feeding expands the region 
of parameter space that allows for a stable host-parasitoid equilibrium. Finally, our results show that 
host-feeding causes inefficiency in the parasitoid population, which yields a higher population of hosts per 
generation. This suggests that host-feeding may have limited long-term impact in terms of suppressing 
host levels for biological control applications. 


1 Introduction 

The host-parasitoid dynamic typically involves a vulnerable period during the year in which hosts are sus¬ 
ceptible to attack by parasitoids. There exists a tendency of synovigenic parasitoids to eat the host without 
laying eggs inside mm- Generally, an adult female parasitoid emerges every year with less eggs than she can 
potentially oviposit in her lifespan. Therefore, a parasitoid will feed on hosts to gain the necessary energy 
and tissue to mature additional eggs mm- However, in doing so, the parasitoid loses a potential host as 
the majority of parasitoids kill the host during host-feeding. This results in a high death rate of the host 
population during the vulnerable period, which can have interesting consequences on later generations mm- 
For further reviews on the host-feeding interaction and its biological implications, we refer the interested 
reader to mmm- 

A traditional approach to describe parasitoid-host dynamics is to use discrete time models, such as the 
Nicholson-Bailey model [9j and others considered more recently [TO] jTTj JT2j [13] [14] . Discrete models can 
monitor the change in population density for distinct points in time, such as each year or each generation, 
which agrees with typical life-cycle themes of insects in temperate climates ns ■ Several authors incorporate 
the effects of host-feeding in continuous models, which are suitable to year-round interactions such as those in 
tropical climates. They conclude that host-feeding can have stabilizing effects mm or no effect on stability 
nans is, while Shea et. al. conclude that egg production delay has a destabilizing effect in host-feeding 
models m- Further reviews on past and recent models are provided by Hassell et. al. in [22 and [23] and 
by Murdoch et. al. in ns • 
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In this work, we use the semi-discrete framework |i23] to develop discrete update functions based on 
a continuous host-feeding dynamic. This modeling approach is more mechanistic as compared to early 
phenomenological models because it incorporates specific host-parasitoid interactions during the vulnerable 
period. We conclude that host-feeding has no effect on stability as compared to previous results by Singh 
et. al |25] ; however, host-feeding combined with density-dependent mortality rates and functional responses 
has a stabilizing effect. The paper is organized as follows: we formulate the model in Section[2] in Section]!] 
we consider the simplest case and show that host-feeding alone does not stabilize the Nicliolson-Bailey model 
[2]; in Sections Q] and [5] we incorporate the effects of density-dependent mortality and quadratic functional 
response, respectively, and consider the stability of the system; we conclude with several discussion topics in 
Section [G] 

2 Model Formulation 

Generally, the discrete-time model used to describe the host-parasitoid dynamics is given by 

H t+1 = F(H u P t ) (1) 

P t+1 =G{H tl P t ), (2) 

where H t and Pt are the adult female host and parasitoid densities, respectively, at the beginning of the 
vulnerable period of year t. where t is an integer. Figured] illustrates the life cycle of the host and parasitoids 
in a year long period. During year t, host pupae mature into H t adults. These H t adults lay a particular 
amount of eggs, which eventually mature into RH t larvae at the beginning of the vulnerable stage. Here, 
R > 1 denotes the number of viable eggs produced by each adult host. The time within the vulnerable stage 
is denoted by r, where r £ [0,T]. Time r = 0 denotes the beginning of the vulnerable stage and r = T 
denotes the end. During this time, adult female parasitoids emerge with egg and seek out host larvae to 
oviposit inside them. The parasitoid egg inside the host becomes a juvenile parasitoid that uses the host as 
a resource. At the end of the vulnerable period, a proportion of the host population is infected. This results 
in an initial parasitoid population the next year as juvenile parasitoids emerge from the host to continue 
their life cycle. A number of hosts escape parasitism, and survive to mature into pupae, which pupate over 
the winter creating the next year’s population of adult hosts. To better understand the discrete update of 
each year, we consider a continuous model to describe the interaction of hosts and parasitoids during the 
vulnerable period. Next, we discuss the approach to model this system of interactions. 

2.1 The Semi-Discrete Framework 

A continuous time model is used to describe the dynamics of the interacting host and parasitoid populations 
during the vulnerable period. The update functions F and G of the discrete model depend on the output 
of the continuous model at the end of the vulnerable period each year. We consider the following chemical 
reaction 


P + L -^4 I + P. (3) 

Here, g{ ■) (units: time -1 parasitoid -1 ) is the attack rate of the parasitoids, which represents the instanta¬ 
neous rate at which the hosts are attacked per parasitoid. This function could potentially be dependent on 
the population of hosts, parasitoids, or infected hosts. In general, we write the continuous model as 


dr - g(-)L{r,t)P(T,t) 

(4) 

dr - 9 (-)L(r,t)P(T,t) 

(5) 

dP(r,t ) Q 
dr 

(6) 
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Figure 1: Life cycle of the host and parasitoids in year t. 


where L(t, t), I(r,t), and P(r,t) denote the concentrations of host larvae, infected larvae, and parasitoids, 
respectively, at time r during the vulnerable period in year t. We subject the system to the following initial 
conditions 


L(0,t) = RH t , J(0,i) = 0, P{0,t) = P t . (7) 

These equations are integrated from r = 0 to r = T. Assuming each parasitized host larvae gives rise to k 
adult parasitoids in the next generation, the update functions in the discrete model are 

H t+1 = F{H u P t ):=L(T,t) (8) 

P t +i = G(H t ,P t ) := kI(T,t). (9) 

Once the system is formulated, investigation into the stability region of Equations m and d2J can be carried 
out using the standard Jury conditions [26J . If we let g = c, i.e. a constant parasitoid attack rate, then we 
obtain the following update 


F(H U P t ) = RH t exp (~cP t T) (10) 

G(H t ,P t ) = kRH t [1 - exp(—cP t T)], (11) 

which is the classic Nicholson-Bailey model. This interaction is unstable and over time both populations 
experience diverging oscillations |15] . Several changes in the form of functional responses, susceptibility 
of attack, and density-dependent mortality can be made to stabilize the system as in mm- Next, we 
illustrate how we incorporate host-feeding into the semi-discrete framework and discuss how to obtain the 
update functions. 

2.2 The Semi-Discrete Framework with Host-Feeding 

The continuous model we consider includes the tendency of the parasitoids to feed on host larvae without 
infecting them with an egg Si. This means the parasitoid population has two phases: without egg and 
with egg. If a parasitoid lays an egg inside the larva, an infected larva is produced and the parasitoid is now 
without an egg. This eggless parasitoid must feed on a host larva before gaining another egg; however, if 
a parasitoid has an egg, it will infect the larvae rather than feed. Implicit to our model is the assumption 
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that once a parasitoid feeds on the host, it immediately gains enough energy to produce an egg, i.e., the 
transitional period from eggless to with egg is instantaneous; however, we do discuss consequences of a delay 
in Section ED This approach is motivated by the work of Shea et. al. m, who consider several parasitoid 
phases with an egg maturation delay. A kinetic diagram below describes the process, 


Pi PL — 

4/ + P 0 

(12) 

Po+P — 

4 Pi. 

(13) 


Here, the function g( ■) (units: time -1 parasitoid -1 ) represents two means by which the parasitoids attack 
hosts. The first reaction is the instantaneous rate at which the hosts are parasitized or infected by a mature 
parasitoid. The second reaction represents the rate at which the hosts are devoured by an eggless parasitoid. 
We assume the two rates are equal as this leads to analytic results. In the most general case, both rates 
may potentially depend on the host, parasitoid, or infected host populations, which yields the following 
continuous model 

dL( £ ^ = - 9 {-)(Po(T,t) + Pi(r,t))L(r,t) (14) 

= g(-)-Pi(' r ,f)T(r,f) (15) 

dP \^ = 9 {-)(Pi(r,t) - P 0 (T,t))L(T,t) (16) 

dPl ^’ ^ = ~ 9 {-)(Pi{r,t) - P 0 (T,t))L(T,t), (17) 

where P(r, t), I(r,t), P 0 (r, t), and Pi(t, f) are the density of host larvae; parasitized host larvae; eggless 
parasitoids; and parasitoids with egg, respectively, at time r during the vulnerable period and in year t. We 
subject the system to the following initial conditions 


L(0,t)=RH u /(0,i) = 0, P o (0,f) = 0, P 1 (0,t)=P t . 


(18) 


This initial condition assumes all parasitoids emerge with an egg; however, we also consider a proportion 
of the population that emerges eggless. We discuss the results of this assumption briefly in the conclusion. 
These equations are integrated from r = 0 to r = T. Assuming each parasitized host larvae gives rise to k 
adult parasitoids in the next generation, the yearly update functions are the same as before in Equations 
m an d ®. We note here that P t = Po(0,t) + Pi(0,f) if there is no density-dependent parasitoid mortality. 
For the remainder of the paper, we suppress the dependence on r and t for the functions P, /, Po, and Pi 
for convenience. 


3 Nicholson-Bailey with Host-Feeding 


In this section, we consider the most simple case. We assume the attack rate g{ ■) in Equations (fT2l) and 
(1131) is constant, i.e. g = c. Our continuous time model during the vulnerable period from Equations (1141) 
G2D becomes 


dL 

dr 

dJ 

dr 

dPo 

dr 

dP 1 

dr 


—c(Pq + P\)L 
cP\L 

c(Pi - P 0 )L 
-c(Pi - P 0 )L. 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 
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In this case, we can solve the system for L and / explicitly © and our discrete update from Equations (HD 
and m can be written as 


Ht+i — RH t e 

D kRH t 
-nt+i = —~— 


—cP+T 


( 1 - 


—cP t T 


) m 
’ 2 





(23) 

(24) 


Here, without loss of generality and for the remainder of the paper, we assume T = 1 and k = 1 since these 
parameters have only a scaling effect on the results. Here, T = 1 corresponds to one vulnerable period, 
which is approximately 90 days. The fixed points to this model satisfy 


H* = 


(3 Ini? 
c(R - 1) ’ 


P* = 


In R 


= 0.7921. 


(25) 


As shown in © this model is unstable for all R. Even if we introduce a transition rate from Pq to Pi via 
an egg maturation delay, the model is still unstable. This suggests that host-feeding alone cannot establish 
stability, which is in agreement to previous models mum. 


4 Density-dependent Host Mortality with Host-Feeding 

It is shown in both phenomenological models ;28] [29] and mechanistic approaches [25] that density-dependent 
host mortality can stabilize discrete-time systems. To investigate stability with host-feeding, we implement a 
density-dependent host mortality in Equations (fill) (fTTl) . The kinetic reactions of the system are depicted 
as 


Pi + L I + P 0 

(26) 

P 0 + L Pi 

(27) 

L 9l ^ - \ Death. 

(28) 


The function g i( •) (units: time -1 ) represents the (potentially density-dependent) host mortality rate due to 
causes other than parasitism and host-feeding. We assume the attack rates are constant, g = c, implying a 
linear functional response. We let the host mortality rate depend on the current amount of host larvae, i.e., 
gi = ci L so that density-dependent effects can act simultaneously with parasitism in the continuous-time 
model. We obtain the following system for r e [0,T], 


^ = -c(Pi + P 0 )L - Cl L 2 

(29) 

dl 

(30) 

Tr= CP ' L 

d f = c(Pi-P 0 ) L 

(31) 

f = - C ^~P^ 

(32) 


subject to the same initial conditions as above. Solving this system explicitly m yields the following discrete 
update system 


H t+ i = 


RH t exp(-cPt) 
f(H t ,P t ) 


Pt+x = ^ln [f(H t ,P t )] + ^ - 1 


(33) 

(34) 
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where 


f(H u P t ) = l + c 1 RH t 


1 - exp(-cP t ) 
cP t 


(35) 


We obtain two nontrivial fixed point solutions. The no-parasitoid (NP) fixed point equilibrium is given by 


Hnp ci R ’ P °' 

The stability analysis in [B] shows that the no-parasitoid equilibrium is stable for 

In R ci 

7 c 


(36) 


We can see from this that sufficiently large values of ci /c, which describes the strength of density-dependent 
mortality versus parasitism, stabilizes the no-parasitoid equilibrium. The second equilibrium point charac¬ 
terizes a presence of both hosts and parasitoids, and is given by 


H* = 


exp (aa) - l \ C P* 

1 _ ex p( 2 l L ) I CliT 


P* = 


In(P) - 


7 = 1.5238. 


(37) 


We consider the region of stable solutions for various values of ci /c. An analysis in [B] shows the equilibrium 
point is stable for 

*• < Si < (38) 

c 7 


where z* satisfies the following equation 

(i?[( 72 * + l) -Ini?] (1 - e~ 7Z ‘) + 1 - eT' 2 *} (l + e~ 2 ^) | 1 

2z*e^ z * (1 - Re~7 z *) h 2 


7 + g l 1 — e 27 ) 


e~ lz = 1 . 


(39) 


Figure [2] shows the stability region as compared to Singh et. al.’s [25j model without host-feeding. We can 
conclude from this analysis that coupling host-feeding with density-dependent host mortality yields a larger 
stability region than that of the same model without host-feeding. Overall, this ensures that the host-feeding 
dynamic has a stabilizing effect. We note that the comparison in Figure [2] is between the host-feeding model 
presented here and the Singh et. al. [25] host-mortality model with parasitic attack rate c/2. In this sense, 
the infection rates of the attacking parasitoids, Pi, of both models are comparable. 


5 Quadratic Functional Response with Host-Feeding 

Previous results show that phenomenological update functions with Type II and Type III functional re¬ 
sponses do not stabilize the Nicholson-Bailey model mm- However, using the semi-discrete framework, 
Singh et. al. showed that a quadratic functional response yields a neutrally stable fixed point with period 
27r/ arctan(VP 2 — 1) in the absence of host-feeding [25]. To investigate the effects of host-feeding, we con¬ 
sider a functional response in the attack rate of Equations m and m so that g = cL. Using Equations 
dUD (HD, our continuous model becomes 


f = -c(P 0 + P l} L 2 

(40) 

dl 9 

(41) 

Tr= CP ' L 

f^c(Pi-P 0 )L 2 

(42) 

^-ctPr-W 2 . 

(43) 
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Stability Region with Host Feeding 



Figure 2: Host-feeding increases the stability region of the Nicholson-Bailey model with density- 
dependent mortality. The stability region specified in (1551) . for the discrete time, host mortality model 
(Equations (1551) and (1551) 1 as a function of the strength of density-dependent mortality to parasitism (ci/c) 
and the number of viable eggs per host, R (black line). The stability region is larger with host-feeding, as 
compared to the density-dependent host mortality model without host-feeding (blue line) [251 . 


subject to the same initial conditions. In[Cl we solve this system explicitly for L and I to obtain the following 
discrete yearly update system 


H t+ i = 


RH t 


Pt +1 = 

The fixed point of the system is 


1 + cRH t P t 

RH t - H t+ 1 + jl - exp 


2(RH t -H t+1 

Pt 




(44) 

(45) 


H* = 



P* 


R- 1 

V7 cR’ 


1.5238. 


(46) 


Stability analysis JC|) shows that this fixed point is stable for all R. Figure [5] illustrates a simulation for 
R = 2 and c = 0.01 with and without host-feeding. We can see that without host-feeding, the populations 
oscillate forever, but with host-feeding, the populations settle to a limiting value. Hence, including host¬ 
feeding with a quadratic functional response in both the attack and infection rate stabilized the originally 
neutrally stable fixed point. It should also be noted that the above host equilibrium is higher than the 
corresponding quadratic functional response equilibrium in [251 with no host-feeding, by a factor of 1.23. 


6 Discussion 

In this paper, we have considered the classic host-parasitoid interaction with a host-feeding dynamic. For 
simplicity and to obtain analytical results, our model ignores gut capacity and includes two states, eggless and 
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Figure 3: Host-feeding stabilizes the Nicholson-Bailey model with quadratic functional response. 

Comparison of trajectories of the dynamic interaction given by Equations (l40l) (l43l) (left) to Singh et. al.’s 

quadratic functional response model without host-feeding [25] (right). Host-feeding with a quadratic func¬ 
tional response stabilizes the oscillatory behavior of Singh et. al.’s model. Simulations are run with R = 2, 
c = 0.01 and initial densities are taken as 8 and 5 for the host and parasitoid population, respectively. 


with egg, where eggs are produced immediately after consumption. In contrast to previous phenomenological 
models, we incorporate the semi-discrete framework, which has more relevance to parasitoid populations with 
one year life cycles and allows us to track the change in hosts during the vulnerable period. The preceding 
analyses show that the effects of host-feeding alone cannot stabilize the classic Nicholson-Bailey model. 
However, including density-dependence and a quadratic functional response coupled with the host-feeding 
dynamic provides a more stabilizing effect as compared to Singh et. al.’s results |25j . For example, when 
R = 2 a stable liost-parasitoid equilibrium exists in a wider range of density-dependent mortality rates by 
a factor of 1.41 so that for higher values of density-dependent mortality, stability still occurs. Indeed, in 
this case, the host-feeding mortality rate is Ci = 0.4549 and without host-feeding, the mortality rate is 
Ci = 0.3466. This means host-feeding relaxes the effect of density-dependent mortality and allows higher 
rates of mortality to exist in the interaction. Furthermore, the mean host density is increased by a factor 
of 1.54 in the host-parasitoid stability region. Therefore, host-feeding makes the parasitoid less efficient in 
reducing the number of hosts in every generation. This may seem contradictory since host-feeding can be 
viewed as an extra process that eliminates the hosts; however, the parasitoid feeds on what could have been a 
viable host for reproduction in the next generation. Therefore, host-feeding may have a short term biological 
control effect but ultimately it reduces the parasitoid population in the long term. In the case of a quadratic 
functional response in the attack rates, Singh et. al. show that the sole equilibrium point is neutrally stable 
[25] . In our model, host-feeding stabilizes this equilibrium point and also shows an increase in the mean host 
density by a factor of 1.23. Hence, host-feeding brings stability to both systems by decreasing the efficiency 
of the parasitoid to reproduce from year to year. 

Investigating other key assumptions in our model may lead to further insight into the host-feeding dy¬ 
namic. For instance, we assume parasitoids emerge into the vulnerable period each year with an egg. As 
Jervis et. al. explain in [32 , the adult parasitoid female emerges with eggs, and later feeding on hosts to gain 
further eggs. However, do all females enter the vulnerable period with eggs? Dieckhoff et. al. and others 






































conclude that many psychological and behavioral factors including resorption and nutrition may contribute 
to the specific egg load at any given time [33]. This translates to the initial condition of our continuous 
model. We investigate this by considering a proportion of the starting population that is initially eggless, 
i.e., Po(0) = AP t and -Pl(O) = (1 — A )P t . We found that as A —> 1 (poor conditions for full egg load at 
adult emergence), both models are less stable. Indeed, for A = 1/2, it can be shown that Equations (121)1) 
(1521) yield an identical update to Singh et. al.’s less stable, host-mortality model, which is the basis for our 
comparison in Figure [2] 

Another investigation we considered is a density-dependent parasitoid mortality. Similar to the model 
presented in Section [I] we consider a parasitoid death rate that is dependent on the total parasitoid popu¬ 
lation. In this sense, our reaction scheme is 


Pi + L 4 / + P 0 

(47) 

P 0 + L 4 P 1 

(48) 

Po - > Death 

(49) 

D C 2 OP 0 +.P 1 ) -p. 

Pi - > Death, 

(50) 


where c represents the constant attack rates by parasitoids as before and C 2 is a constant. This model is 
analogous to Singh et. al.’s density-dependent, parasitoid mortality model without host-feeding [25]. As 
shown by Singh et. ah, the non-host-feeding system is stable if C 2 > c. In the host-feeding case, an analytical 
solution for I(r,t) cannot be obtained; however, numerical investigation suggests that the system is stable 
for C 2 values that are slightly less than c. Indeed, a trajectory with R = 2, c = 0.01, and C 2 = 0.0098 
yields asymptotically stable results whereas the model without host-feeding experiences diverging oscillations. 
Although the stability region is expanded modestly, we can ultimately conclude that host-feeding has a 
stabilizing effect in the density-dependent parasitoid mortality case. 

The study by Shea et. al. [2Tj focused on the effects of egg limitation [34] in a host-feeding interaction. 
They conclude that stability is effected by the length of latent period, i.e., the time it takes for eggs to 
mature as the gut is emptied. Indeed, the longer the latent period, the longer it takes for the system to 
stabilize. Among other conclusions, they also observe that stability is not effected by the number of eggs that 
can be stored. In this sense, our model considers an instantaneous egg maturation process and parasitoids 
carry a single egg. Because we only consider a single egg, we are able to study the system analytically and 
gather the results shown above. However, to discuss the case of egg maturation, we implement a delay in 
egg production into the semi-discrete framework with the following reactions 


Pi+L 

^I + Po 

(51) 

Po+L 

Pl/2 

(52) 

Pi ° 3 

Pl/2 - 

► Pi. 

(53) 


As compared to the original model in Section [2J we add a population, P 1 / 2 , that is analogous to Shea 
et. al.’s P 01 population, which represents a parasitoid with zero eggs and a full gut. The constant rate C 3 
measures the transition from eggless to with egg as the parasitoid gains enough energy to effectively mature 
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its egg. Using these reactions, our equations become 


^ = -c(P 0 + Pi)L 2 
dr 

(54) 

T = cPlL2 

dr 

(55) 

- w 2 

(56) 

dPi/o 0 

= cPiL 2 - C 3 P 1/2 

(57) 

= C 3 P 1/2 ~ cP\L 2 . 

(58) 


Analytical solutions to this model are unattainable, but we confirm the destabilizing effect of the latent 
period numerically. In Figure U we see a series of host densities using the quadratic functional response 
model plotted against time with a varying delay parameter, C 3 . As the rate of maturation gets slower, the 
delay increases. For a critical delay the stable host-feeding, functional response model becomes unstable as 
the density of hosts begins to experience diverging oscillations. 

It is worthwhile to note that the characteristic time scale, T, in which the vulnerable period takes place 
is approximately 90 days. The instability is caused by a delay that is approximately 0.18, which suggests 
that the egg maturation delay is approximately 16 days, which is far too long for an egg to mature. As 
illustrated in Figure [I] a delay value closer to 10~ 2 or 10~ 3 , which means it takes on the order of hours for 
an egg to mature, yields stable results. Therefore, in this more realistic case, the system is always stable. 


Delayed Functional Response 



Figure 4: Incorporating a delay due to egg production has a destabilizing effect on the quadratic 
functional response model. A series of host population density ( 2 -axis) trajectories are plotted against 
time (x-axis) for each value of egg maturation delay (y-axis). As the delay time, I/C 3 , approaches 0.18 or 
approximately 16 days, the system becomes unstable. However, the system is always stable for smaller (on 
the order of hours) delay time. Simulations are run using continuous system given by Equations (l54|) - (l58l) 
and discrete yearly update given by Equations © and © with T = 1 and k = 1. Parameters: R = 2, 
c = 0.01, H( 0) = 8, and P( 0) = 5. 
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Finally, we note that many other important consumer-resource dynamics can be implemented into the 
semi-discrete framework. A necessary future step would be to investigate the effects of a larger egg load 
and/or to implement a probability that a parasitoid will host feed. In this sense, an optimal stability 
criterion could exist. Furthermore, we seek to investigate susceptibility of host risk to parasitoid attack, as 
in m- Susceptibility of risk cannot stabilize the model with host-feeding alone, but could provide interesting 
results with host-dependent mortality or a quadratic functional response. The generality of the semi-discrete 
framework allows us to implement these changes with ease. 
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dPp dPi _ 
dr dr 


A Analysis of Density Independent Mortality 

We consider the explicit solution to Equations © m- Adding ED and (l22l) gives 

P 0 + Pi = P t . (59) 

Substituting Pi = P t — P 0 into m yields 

=>- L(r,t) = RH t ex p (— cPtT ). (60) 

We can now solve for P 0 by substituting our expression for L and Pi into equation ED to obtain 

(61) 

Using the integrating factor to solve this equation, we get an expression for both P 0 and Pi. We have 

(62) 


dL 

— = —cP t L 
dr 


dP 0 


— -h 2 cRH t exp(-cP t r) P 0 = cRH t P t ex p (~cP t r). 

ar 


n / A R t Pt 

p o(r,i) = Y - Y 


Pi 


exp | [l - exp (-cP t r) ] j 

(r,t) = y + y exp j [l - exp {-cP t T) ] j . 


(63) 


(64) 


Finally, substituting our expression for Pi into equation (1201) gives 

dl cRH t P t , n , cRH t P t r r, \ f RHf ri t D ill 

— = -^-exp (-ci? t r) +- - -exp (-cP t r) exp j —- [1 - exp (-cP t r) J | . 

Solving this and applying the initial condition gives the following expression for I (r, t) 

[l ^ ex P (—cPtr) ] + y jexp j^-[l - exp (-cP 4 t) ] | - l|. (65) 

Using the solutions for L and /, we find the update as in Equations (B^l) and (HH1) using the definition in 
Equations © and ©■ We can analyze the resulting discrete-time model by implementing the procedure 
outlined in (26]. Using a general discrete model such as 


H t+1 =F(H t ,P t ) 
Pt +i = G(H t ,Pt), 


( 66 ) 

(67) 


we can perform a linear stability analysis about the fixed point (i7*,P*), where H* and P* satisfy the 
following system of equations 


H* = F(H* 7 P*) 
P* = G(H*,P*), 


( 68 ) 

(69) 


by determining if the spectral radius of the Jacobian matrix evaluated at the fixed point is less than one. 
That is, if the magnitude of the eigenvalues of the following Jacobian matrix, 


J := J(H*,P*) = 


" dF 


dF 


dH 

(■ H*,P *) 

dP 

(H* ,P*) 

dG 

dH 


dG 

dP 



(70) 


are within the unit circle, then the fixed point (H*,P*) is asymptotically stable. The eigenvalues fall within 
the unit circle if the following three Jury conditions hold, 


1 — Tr(J) + Det(J) > 0 

(71) 

1 + Tr(J) + Det(J) > 0 

(72) 

1 - Det(J) > 0. 

(73) 
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In the case of density independent mortality, our update functions are given by Equations El) and (El) . 


F(H t ,P t ) = RH t e~ cPt 
G(H t ,P t ) = ^( l-e- cP ‘) 


P 

2 


/ 

\RH t , 


< exp 

p, ( 

1 — e f ) 



with fixed point 


pin R 
c (R-iy 

Using these expressions, we can evaluate the trace 


Tr (J) 
Det(J) 


P{l + eP) 

2 

(3(1 + e 0 ) 
2 


P* = —, p = 0.7921. 
c 

and determinant of the Jacobian matrix 

e^(l — P) + 1 Ini? 

P(1 + e@) + R- 1 

"e^(l — /?) — 1 , Ini? ' 

/?(1 + e@) + R- 1 ' 


as 


(74) 

(75) 

(76) 

(77) 

(78) 


We can see that 1 — Tr(J) + Det(J) = 0 for all R. 


Hence, the fixed point is not asymptotically stable. 


B Analysis of density-dependent Mortality Model 

We consider the explicit solution to Equations (l29l) (l32l) . Adding (IdTl) and (l32l) gives Pq + P\ = Pt, which 

upon substitution of Pi = Pt — Po into Equation (l2fll) yields 


dL 

dr 


(—cPt + C\L) L 


L(t, t) 


RH t exp(-cP t r) 

1 + Cl RH t 1 - exp e ( ~ cPtr) ' 


(79) 


The previous differential equation can be solved using partial fraction decomposition. We can now solve for 
Po by substituting our expression for L and Pl into equation (l3lTl to obtain 


dP 0 2cRH t exp (-cP t r) ^ _ cRH t P t ex p (~cP t r) 
dr i + Cl RH t 1 ~ exp e ( ^ cPtT) ° _ 1 + CiRH t 1 ~ exp e ^ cPtT) ' 


(80) 


Using the integrating factor to solve this equation, we get an expression for both 
explicitly for Pp and using the initial condition from the previous section yields our 
Pi as 


Po{r,t) 


P 

2 


1 


2 


1 + ciRH t 


1— exp( — cP t r) 

cPt 


2 c 
C 1 


Po and Pl- Solving this 
solution for both Po and 

(81) 


Pi(r,f) 


P 

2 


1 + c\RHt 


1—exp(— cP t r) 
cPt 


Finally, substituting our expression for Pi into equation (1801) gives 


(82) 


dJ 

dr 


2 


cP t exp (~cP t r) 

1 + ciRHt 1 - exp e ( ~ cPtT) 


cexp(— cPtr) 


1 + ciRHt 1 - exp e ^ cPtT) 


Hi 


(83) 


Solving this and applying the initial condition gives the following expression for I (r, t) 




1 + ciRH t 


1 - exp(—cPtr) 


cPt 


P 

4 


1 + CiRHt 


1 - exp(—cP t r) 

cPt 


- 1 


(84) 
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Using the solutions for L and /, we find the update as in Equations (1551) and (1551) using the definition in 
Equations (0 and 0. In the density-dependent host mortality case, our discrete update functions are, 


F(H t , P t ) = 


RH t exp(— cP t ) 


G(H t ,P t ) = ^ In 


1 + c± RHt 


1 - exp(-cP t ) 

cPt 


P 

4 


1 + c\RH t 


1 - exp (~cP t ) 
cPt 


-n. 


(85) 


( 86 ) 


This system has two nontrivial fixed points. As noted in the main text, the first no-parasitoid equilibrium 
point is given by 

h ng = P * = 0. (87) 

Denoting the Jacobian matrix for the first fixed point by Jnc-, we find the trace and determinant to be 

2 In P 


Tr( Ao , = I-I 


Det(Jjvc) = — 


1 

4P 


R~ - 1 - 
R~i - 1 - 


z 

2 In R 


( 88 ) 

(89) 


where z = C\/c. Using these expressions, Jury condition 2 gives the most strict condition so that z > In P /7 
for stability. The second equilibrium point is given by 


H* = 


exp iff) - 1 \ C P* 


1 - 


exp(2gl) 


P* = 


In(P) - 


Cl 


7 = 1.5238. 


(90) 


Define J as the Jacobian matrix evaluated at the second equilibrium point above, then the trace and deter¬ 
minant are given by 


Tr(J) = 
Det(J) = 


{e^ z [In R — (jz - 1 )] -R}{ 1 - e l ' 2 )(l + e~ 2 f 1 


2 z(R-e'i z ) 2 

{R [In R - ( 7 * + 1)] (1 - e~ lz ) + - 1} (l + e~ 2 f 1 

2 z(R-e^ z ) + 2 


7+i(l-e-7 




7 + ^ (l — e -27 ) 


o~~i z 


(91) 

(92) 


Here, Jury condition 3 is the most strict, and we must have z* < z < In R/ 7 , where z* solves 1 — Det( J) = 0. 
Figure [5] shows a plot for both stability regions. 


C Analysis of Quadratic Functional Response Model 

We consider the explicit solution to Equations (1401) (1431) . Adding (1401) and (155|) gives P 0 + Pi = P t - 

Substituting Pi into Equation (1501) yields 


dL 9 

— = -cP t L 2 
dr 


L(t, t) = 


RHt 


1 + cRH t PtT 

We can now solve for P 0 by substituting our expression for L and Pi into equation (1501) to obtain 


(93) 


rP, = 


H r t 


dr ' (l + cPP t P t T) 2 ^ u (1 + cRH t P t T) 2 ' 

Using the integrating factor to solve this equation, we get an expression for both P 0 and Pi. We have 

'2 RH t ( 1 


p / ,N Pt Pt 
Po{r,t) = — - — exp 


1 + cRH t PtJ 


- 1 


(94) 


(95) 
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p / Pt Pt 

Pi{T,t) = — + — exp 


2 RH t 


1 


1 + cRH t P t T 


- 1 


Substituting our expressions for Pi and L into Equation (14T1) . we obtain the following ODE for /, 

cR 2 H 2 P t exp (-2^) exp f 2RHt 


dl 


cR 2 H 2 P t 


P t (l+cRH t P t r) 


dr 2(1 + cRH t Ptr) 2 2(1 + cPP t P t r) 2 

We can solve this equation directly to obtain the following expression for I, 


I(r,t) = 


RH t 


1 - 


1 


1 + cRH t P t r 


Pt 

— < 1 — exp 


2 RH t 


1 


Pt V1 + cRH t Ptr 


- 1 


(96) 


(97) 


(98) 


Using the solutions for L and /, we find the update as in Equations (l44l) and (l45l) using the definition in 
Equations (0 and 0. In the quadratic functional response case, our discrete update functions are, 


F(H t ,P t ) = 
G(H t ,P t ) = 


RH t 


1 + cRH t P t 
RH t 


1 - 


1 


2 V 1 + cRH t Pt 

which gives rise to the following fixed point 


4 


1 — exp 


2 RH t 
Pt 


( --- 

\1 +cRHtPt 



H* 

The trace and determinant of the 

Tr(J) 

Det(J) 

All three Jury conditions hold for 


[T p ._ i 

V cR 1 y/j cP’ 


Jacobian matrix is given by 



= i ( 1 + e-*) + 5 L[l -( 1 + 


R > 1. 


7 = 1.5238. 

-|( 1 + e- 27 ) 

27)e -27 ] ■ 


(99) 

( 100 ) 

( 101 ) 

( 102 ) 

(103) 
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